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^ . It is now a known fact that if we happen to be living in the middle of a large underdense region, 

then we will observe an "apparent acceleration", even when any form of dark energy is absent. In 
I this paper, we present a "Minimal Void" scenario, i.e. a "void" with minimal underdensity contrast 

(of about —0.4) and radius (~ 200 — 250 Mpc//i) that can, not only explain the supernovae data, 
but also be consistent with the 3-yr WMAP data. We also discuss consistency of our model with 
various other measurements such as Big Bang Nucleosynthesis, Baryon Acoustic Oscillations and 
, local measurements of the Hubble parameter, and also point out possible observable signatures. 

o 
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■ One of the most baffling problems in cosmology and fundamental physics today concerns the acceleration of the 
, universe, as inferred from the redshifting of the type la supernovae. Along with this observation, the WMAP data 
• and the large scale structure measurements can all be explained by invoking a dark fluid with negative pressure dubbed 
\ as dark energy. This has given rise to the so-called flat ACDM or concordance model consisting of approximately only 

4% of visible matter (baryons), the rest being dark (approximately 3/4 dark energy and 1/4 dark matter). However, 

J — what is this dark energy and why its abundance should be such that it happens to be exactly in concordance with 

, matter density today, remains very much a mystery. 

J> . Recently, a few researchers have tried to take a different point of view: what if the effect of large scale structure 
could account for the observed luminosity to redshift behavior of type la supernovae {i.e. give rise to an "apparent" 
acceleration of the universe), without Dark Energy? This question is important because an affirmative answer might 

■ obviate the need for a dark energy component/cosmological constant, which has presented a plethora of unresolved 
theoretical issues. Recent studies of exact solutions to the Einstein equations have, in fact, been able to reproduce the 
observed luminosity to redshift relation that is usually attributed to acceleration, provided that we lived in a large 
region ( "void" ) that has less matter density than the spatial average density over the cosmological Horizon [l], 0, Q 
(see for a review). One might naively conclude that this result can obviate the need for dark energy. However, in 
order for the void model to be taken seriously, several key issues have to be addressed. 

Firstly, the observation of small, nearly scale invariant CMB temperature fluctuations, strongly supports the principle 
that our universe is homogeneous and isotropic on large scales. In our present Universe non-linear large scale structures 
exist, marking a deviation from homogeneity; however, according to our current understanding of structure formation, 
0(1) non-linearities are only expected typically at scales ~ O(10Mpc//i). In this case one can again argue that the 



I. INTRODUCTION 



'Electronic address: stephon@slac . Stanford . edu | 
t Electronic address: 'tbi^was^gravitylpsu'edi? 
-^Electronic address: notari@hep. physics. mcgill .ca] 
§ Electronic address: |deepak@pliy s . psu . edu | 



effect of these inhomogeneities on cosmology, which is governed by the Hubble scale ^ 3000Mpc//i, would be to too 
small to be significant. However, there are reasons why one could be wary of such a conclusion. 

From the theoretical point of view, the non-linear behaviour of structure formation is not a trivial issue. For instance, 
due to non-linear effects it is known that smaller voids can percolate to form much larger underdense structures which 
occupy most of the volume of the Universe (see e.g. according to which such a percolation has a threshold, when 
the density is about 50% lower than the average), forming what is known as a "cosmic web" of superclusters and 
voids. Also, we note that non-standard features on the primordial power spectrum, such as a spike at a particular 
scale, or some non gaussianity may enhance the possibility of having larger structures and voids. 

Observationally speaking, several huge nonlinear structures (notably, the Sloan Great Wall has a length of 400//i 
Mpc %) have been revealed through surveys like SDSS and 2dF (of course, these data are only tracing the visible 
matter, so their interpretation in terms of total matter is subject to a bias). It is unclear whether the presence of these 
large observed objects is consistent with the present understanding of structure formation. For example Einasto Q 
claims a discrepancy (by a factor of 5) between the observed abundance of such objects (superclusters) and the values 
obtained using N-body numerical simulations. Peebles has also argued that our understanding of structure formation 
and observed voids are in apparent contradiction Q, and that this can be classified as a 'crisis' of the ACDM model. 
Further, there has been observational evidence for the presence of a local large underdense (~ 25% less dense) region 
(that extends to ^ 200Mpc/ft,) from number counts of galaxies Q. This represents a 4 sigma fluctuation, and would 
be at odds with ACDM. More recently, there has been a claim that the presence of the cold spot in the CMB detected 
in the WMAP sky [13] is also associated with a similar Big Void in the large scale-structure Intriguingly, the 
presence of such Big Voids has also been advocated by in order to explain some features of the low multipole 
anomalies in the CMB data (in addition to the cold spot). Finally we note that two recent papers [H, Ql claim a 
significant (95% C.LO detection of an anisotropy in the local Hubble flow in the Hubble Key Project data [lj| and in 
the SN la dataset [13|. This would be a completely natural consequence of being inside a large local void since, 
of course, we are not expected to be exactly at the center and the void is not expected to be exactly spherical. 

To summarize, the large scale structure of our universe might be richer than we thought, which can have far reaching 
consequences for cosmology. However, it is fair to say that the presence of large voids becomes more unlikely (thus 
requiring probably a non-conventional paradigm of structure formation), as the size of the void and the density 
contrast that we consider become larger. This emphasizes the need to find the "Minimal Void (MV) Model" i.e., 
with minimal length scale and underdensity contrast that is required to give a consistent fit to the supernovae data 
(the reader will easily recognize that the larger is the void, in general the better is the fit) . This is the first goal of 
our paper. We find that to obtain an acceptable fit (with goodness-of-fit^ close to 50%) to the current supernovae 
data one needs "us" to be located roughly centrally (with 10% precision in the radial position) within an underdense 
region stretching upto a redshift ~ 0.08. If one is willing to live with a worse fit (goodness-of-fit ^ 10%) then one 
can even go down to z ~ 0.055. The underdensity that is needed is of about 5 ~ —0.4. Now, this is a very large 
region (corresponding to a radius between I6O//1 Mpc and 250/ h Mpc). However we believe that finding a viable 
alternative to the presence of Dark Energy is a task which is important enough to consider such possibilities (and, as 
noted before, similar structures have been advocated for solving other problems in cosmology, as the low-Z anomalies 
and the cold spot in the CMB). 

As an aside we note that we obtain analytical expressions for the luminosity-redshift curve for arbitrary density profiles 
which are excellent approximations even when the local inhomogeneous patch extends up to ~ 400 Mpc/h. 

The second important issue that one has to address in the context of the MV model is whether it can reproduce 
the successes of ACDM model for many different observations, most notably the WMAP third year data. In this 
paper we present an analysis of the MV model subject to the WMAP Syr data using the COSMOMC package 
developed to perform a likelihood analysis of theoretical parameters using a Monte Carlo Markov Chain (MCMC) 
method; we refine the analysis of the type la supernovae data and we combine them together. We find that using 
standard statistical analysis, the MV model accounts for the WMAP and SN la data while being consistent also with 
local measurements of the Hubble parameter. 

In what follows, we will find a consistent fit to both the WMAP and SN la data, provided the Hubble parameter 
ifout = /iout/3000Mpc~^ outside the void is very low, hout ^ 0.45. Then the void plays the role of providing a higher 
value for the local measurements of the Hubble parameter (iJin = /i/3000 Mpc~^). It is exactly this jump in the 



^ The goodness-of-fit for a fit is tlie probability that, given a set of physical parameters, the data has a smaller or equal than the 
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Hubble parameter that gives rise to an apparent acceleration. Additionally, in order to fit the CMB, the primordial 
spectrum has to deviate from the usual nearly flat spectrum. Specifically we try the fit allowing for running of the 
spectral index in the observed 7 e-folds of the CMB sky. Our best-fit has a low spectral index with a large running. 
The overall goodness-of-fit to the WMAP 3 yr data for our best-fit model is around 26% as compared to 41% of the 
ACDM model. 

We should clarify that although we quote comparative statistics between MV and ACDM model, it is only meant as a 
guide, our aim here is not to compete with the ACDM model. According to the Bayesian statistical likelihood analysis 
of both the supernovae and the CMB data, our best fit MV model is still disfavored by many standard deviations 
as compared to the concordant ACDM model. Crucially however, such an inference is based on assuming a "flat" 
prior on the value of the cosmological constant. In other words it relies on the a priori assumption that all the values 
of the cosmological constant are equally likely. According to the Bayesian theory, such a priori probablities are to 
be assigned based on theoretical prejudice. Unfortunately our understanding of the cosmological constant is rather 
incomplete to say the least. As is well known, theoretical expectations suggest an enormously large value ~ Af^, and 
even with supersymmetry it's "natural" value should have been around (TeV)*, in obvious disagreement with our 
universe. Accordingly, before the discovery of our accelerated expansion, our theoretical prejudice had been to assume 
that the cosmological constant must in fact vanish possibly due to some symmetry or other theoretical considerations 
(for a recent review see for instance ^20]). Here we take the same approach, that the "flat prior" assumption may 
actually be misleading and therefore a direct likelihood comparison between a A ^ model with a A = model may 
not be appropriate. Rather we should focus on "independent" statistical quantities such as "goodness of fit" which 
can simply test whether a given theoretical model is consistent with the observational data. In other words, if we had 
a different theoretical prejudice (for example that a non-zero cosmological constant is "unphysical" ) , then we could 
just ask the question whether a non-homogenous matter distribution can fit the data, with an acceptable value of the 
goodness-of-fit. To summarize, although our model has a worse fit than ACDM, in our opinion the statistics suggest 
that our void model in an EdS background can still be consistent with SN and CMB. 

It is natural though to wonder whether one can make these fits better by including perhaps more parameters. We 
consider two such possibilities in brief. In the authors obtained a slightly better fit to the WMAP data, without 
Dark Energy, as compared to the ACDM model by including a bump in the spectrum at some particular scale (see 
also [H, [I^l). We discuss how this can be integrated in the MV framework. Moreover, this idea looks particularly 
appealing, because the existence of a bump in the primordial spectrum could, in fact, enhance the probability of 
finding large voids in the present Universe (the scale of the bump happens to be roughly the scale that we need for 
a Minimal Void). Next, we consider the possibility of having a slight curvature in the model. It turns out that this 
also improves the fit to WMAP considerably. 

Let us now come to the question of consistency between our model and the measurements of the local Hubble 
parameter. Although different observations suggest rather different values, 

0.55 </i< 0.8, (1) 

is perhaps a fair range to consider. As we will find out, the supernovae data essentially constrains the amount of 
jump, hont/h (or equivalently the underdensity contrast in the void), to some range. Combining this with the WMAP 
analysis (which constrains /lout) we get an allowed range for h. These allowed values are definitely low, but we find 
that our h can be as high as 0.59 and therefore be consistent with the local measurements of the Hubble parameter. 

Finally, we briefly discuss consistency of our model with various other measurements, such as baryon acoustic os- 
cillations (BAO), baryon density obtained from Big Bang Nucleosynthesis (BBN), constraints on as coming from 
weak lensing experiments. Integrated Sachs- Wolfe (ISW) effect, etc. An important task that we leave for future is to 
perform an analysis of the SDSS data including Lyman-a forests, without which one cannot really pronounce the MV 
model as a viable alternative to the concordant ACDM model. 

We now proceed as follows: in section [Hi we introduce our swiss-cheese model and briefly discuss the non-linear 
structure formation captured in this model, as well as the photon propagation in this configuration. In section IIIII 
we explain qualitatively how the MV model can be consistent with both the supernovae and WMAP data, as well 
as local measurements of the Hubble parameter. In section IIVI we perform supernovae fits for the void model. 
This includes finding a SN-I best-fit parameter set and comparing it with values for the ACDM model, as well 
as finding a combined best-fit parameter set, which has the maximal jump (this will be needed to better fit the 
WMAP data) with "acceptable" ■ Next in section |Vl we perform a MCMC analysis of the WMAP data without a 
cosmological constant. Again, this involves obtaining a WMAP best-fit parameter set, and also finding a Combined 
best-fit model consistent with supernovae with reasonable as compared to the best fit "concordance" ACDM model. 
In section fVlH we briefly discuss consistency of MV model with other observations such as BBN and BAO. Finally, 
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we conclude summarizing our findings and also pointing out unique predictions of the MV model. The appendix 
contains approximate analytical solution of the trajectory, redshift and luminosity distance of a photon in the radially 
inhomogeneous "LTB" (Lemaitre-Tolman-Biondi) metric. 



II. LARGE SCALE STRUCTURE AND LTB METRICS 

As emphasized in the introduction, we are currently living in a universe with significant inhomogcneitics: non-linear 
structures and voids are expected on average at scales ~ 0(10) Mpc//i, and there is observation of structures up to 
much larger scales, ~ 300Mpc//i. In this paper we will advocate that perhaps we are sitting near the centre of a "Big 
Void" spanning a radius of ~ 200 Mpc//i which, as we will explain, is roughly the minimal size needed to account for 
the SN-Ia supernova data (although one can go down to values of about 150Mpc//i by accepting a slightly worse fit). 

An accurate way to model such inhomogeneous structures/ voids, which avoids any possible pit-falls of perturbative 
arguments, is to use exact solutions of General Relativity that can be studied both analytically and numerically. In 
particular we focus on spherically symmetric LTB metrics [2l| to describe "radially" inhomogeneous patches of any 
desired radius, L (such metric describes the most generic spherically symmetric dust-filled spacetime; we refer to 
appendix IIXI for definitions and details). Such spherical patches can be pasted onto a homogeneous FLRW metric 
consistently [22]. It also ensures that the average density inside the spherical patch is the same (almost exactly, see 
again appendix IIXI for details) as the background density outside the patch. Thus an underdensity around the central 
region is compensated by a shell-like structure near the circumference^. 

Technically, it is somewhat complicated to describe the dynamics of the LTB metric (see appendix IIXI for details and 
for the choice we made for the so-called mass function), but intuitively it is as if one had an independent scale factor 
corresponding to each (comoving) radial coordinate, r, which is evolving as an independent FLRW metric with a given 
spatial curvature k{r). A priori, k{r) is an arbitrary function which also determines the density profile. Assuming 
L <C Rh (the Hubble radius) one has 

(p) (t) M2 

^(^'^)" l + (,/t)^/3,(,) - -here (,)(t)^^, and e(r) ^ 3fc(r) -I- rfc'(r) . (2) 

We observe that the FLRW behaviour for the density is given by the factor (p) (t) , while the fluctuations are provided by 
the presence of e(r) in the denominator. When e(r) is close to its maximum value we have a void, while when it is close 
to its minimum, it signals an overdensity. Note that at early times the density contrast (S(r, t) = {p{r, t)— (p) {t))/{p){t), 
defined in the usual way, grows as t'^^^, in agreement with the prediction of cosmological perturbation theory. On the 
other hand at late times, when {t/tQ)^^^e{r) ~ C(l), the density contrast grows rapidly (and this result is the same as 
found within the Zeldovich approximation [23]). In fact, for an overdense region, the structure ultimately collapses, 
as to be expected because LTB metrics cannot account for virialization that we observe in nature. Nevertheless, for 
our purpose, as long as we do not reach the collapse time, LTB metrics adequately capture the effects of non-linear 
stucture formation on photon propagation. 

Now, we are interested in modeling a spherical void region surrounded by a compensating shell-like structure, and 
this is obtained using a k{r) which starts off from a maximum at r = and falls off to a constant value at r = L such 
that 

k'{Q) = k'{L) = 0, (3) 
k{L) = ^Uk, for]r!fc]«l, (4) 

One can check that such an LTB metric can consistently match to an FLRW background [13], with curvature abun- 
dance Q.k- In this paper we will mostly focus on a background FLRW metric which is flat. The essential reason for 
choosing a flat background metric is that curvature is known to be constrained to be very small in order to get a good 
fit of the WMAP data along with other measurements (such as measurements of the Hubble constant [11]). However, 
in section ?? we will present a brief discussion on how things may change in the context on the MV model if we allow 



^ In fact we may speculate that the Great Sloan Wall may be indicative of such a shell-like structure, given its location, at about 250 Mpc/h 
away from us, and its two-dimensional shape [^. 



4 



for curvature. We note in passing that in LTB models we are considering we do not have back-reaction effects in the 
outside region, i.e. on the average the FLRW regions do not feel at all the presence of holes. The particular choice of 
the curvature function that we employ to model the inhomogeneities and fit the supernova data is given by 



k{r) = fc„ 



(5) 



One can check that Eq.Q satisfies Eq.(j4]), in the case with ^lk = 0. It contains two important physical parameters, L 
and /Emax, which correspond to the length-scale and amplitude of fluctuations respectively ^. In the rest of the paper, 
this is the profile that we will focus on, although some of the analytical results are general for any fc(r). 



III. THE MINIMAL VOID MODEL 



By now in a scries of papers [l|, 0, Q it has been shown that a "large" local underdcnsity can reproduce reasonably 
well the luminosity distance versus redshift, Dl(z), curve that one observes, and thereby can mimic dark energy (for 
slightly different approaches based on inhomogeneities, see [13, HE, |26„ J8J ) . However, the reason one is skeptical 
of such an explanation is because a straightforward extrapolation of the density fluctuations observed in CMB gives 
us today a scale on nonlinearity (that is, the scale in which the expected density contrast is of 0(1)) of at most 
^ O{10)/h Mpc, much too small to explain away dark energy; as we shall see later, we need to invoke a Big Void with 
a radius of about 200//i Mpc (and with average density contrast of roug hly ((5^) ~ 0.4). The probabihty of having 
non-linear structures at larger scales becomes progressively smaller. Using the conventional linear and Gaussian power 
spectrum for radii of about ~ 200 Mpc/h the typical density contrast instead is only of about 0.03 — 0.05 (for a radius 
of ~ 160 Mpc/ft, the typical contrast is instead about 0.06). However, as argued in the introduction, one cannot take 
such an analysis at its face value. There are both theoretical and observational suggestions that we might actually 
have larger underdensities in such voids in our universe. 

Nevertheless, it is clear that the presence of large voids becomes more and more unlikely (or that it would require a 
non-conventional paradigm of structure formation) as the size of the void and the density contrast become larger and 
larger. 

This emphasizes the need to find the "Minimal Void Model" i.e., with minimal length scale and underdcnsity contrast 
that is required to give a consistent fit to the supernova data. This is obtained by realizing that the crucial evidence 
for acceleration comes from the fact that we observe a mismatch between the expansion at low redshifts (between 
roughly 0.03 < z < 0.07) and the expansion at higher redshifts (where supernovae are observed [1^, between roughly 
0.4 < z < 1). This situation arises because of the current experimental status of supernovae observations: we have 
very few data in the redshift ran ge b etween 0.07 and 0.4 (the situation will dramatically change with the coming release 
of the SDSS-II supernovae data [30|). Thus it is not necessary to alter the EdS Dl{z) all the way up to z ^ C'(l), but 
a large correction to the Hubble expansion in the local region, 0.03 < z < 0.07, stretching up to ~ 200Mpc//i, may be 
sufhcient. In particular if we are living in a local underdcnsity, then we experience extra stretching as voids become 
"more void" (that is how structure formation works) which manifests as a local Hubble expansion rate larger than 
average (outside the patch) , precisely what is required to mimic acceleration. Another way of seeing this is that all 
sources in the local region have a collective radial peculiar velocity due to the gravitational attraction of the shell-like 
structure, which adds to the overall expansion. 

We may also note that recently [si*] has claimed a possible detection of a jump in the supernova Hubble diagram, 
exactly in the direction of having a large void. However the void radius (about 75Mpc/ft,) and the jump (about 7%) 
are smaller than what we are proposing. 

Let us now see more precisely how the MV model works. Let us start with the observation that the Dl{z) corre- 
sponding to ACDM model is in good agreement with the observed supernovae Thus, if we can ensure that our 
MV model can approximately agree with the ACDM Dl{z) curve both in the low and high redshift supernovae range, 
then we can expect to find a good fit to the data as well. We first focus on the high redshift region, i.e. outside the 



^ The exponent of r/L has been chosen to be equal to 4, but the reader may note that any exponent n > 1 is good, as well. Varying n 
one varies the width of the shell-like structure. The larger the n, the flatter the void, and narrower the structure. However, we choose 
to stick only to the case n = 4, since it already gives us a sufficiently flat proflle for the underdense region which we found to improve 
the supernova flt, and anyhow the whole analysis and discussion is not very much affected by the precise shape of the shell. 
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LTB patch. In this region the Dl{z) curve of the MV model basically corresponds to that of the homogeneous EdS 
curve parameterized by the lower average Hubble parameter'', hout- Further, in this range of high redshift supernovae, 
the EdS curve can run very close to the ACDM model, albeit with a different, slightly lower, Hubble parameter, /lout 
as compared to the Hubble parameter h of the ACDM curve. For instance, if we compare the EdS distance (De) 
with the ACDM distance (Da) 

(6) 

it turns out that the ratio TZ does not change much in the relevant range of high-z supernovae, 0.4 < z < I: 

7^(0.4)/7^(l) ~ 1.12. (7) 

Moreover, the ratio TZ{z) itself is proportional to the ratio /i//iout- Thus, by choosing the latter ratio appropriately, 
the luminosity distance of the average EdS model can be made to approximately coincide with that of the ACDM 
one in the redshift range 0.4 < z < 1, and consequently one expects that the EdS/MV model will be consistent with 
the high redshift supernovae. 

Next, let us look at the low redshift region. In this region, the Dl{z) curve is basically linear, the slope being given 
by the Hubble parameter: 

i/-.lim^ .M!£. (8) 

2-»o z h 

Thus in order for the MV model to agree with the best-fit ACDM, the Hubble parameter inside the LTB patch should 
coincide with the measured local Hubble parameter. In other words, if the MV model can account for the jump, J, 
between the locally measured Hubble parameter h inside the patch, and the lower average Hubble parameter, /lout, 
outside the patch: 

^-T^, (9) 

"out 

then we expect to have a good agreement with the supernovae data. 
Thus the challenges are 

• to quantitatively verify our above hypothesis of being able to find a good fit to the supernovae with an appropriate 
"jump" . 

• to find whether local inhomogeneities in an LTB model can account for such jumps. 

Provided we can make this work, such an analysis will also tell us what a good range for the jump parameter is. 

As we will see later, we find a very good fit to the SN data (where we use the dataset '29']), with goodness-of-fit 
~ 50%, without A in the MV models. Assuming our model, a parameter estimation (with likelihood /^) gives at 
95% C.L. the following range for the jump parameter: 

1.17 < < 1.25. (10) 

On the other hand the fit to the WMAP data will fix the value of the global hout- As we will see in section |Vl this is 
the important quantity for the photons that come from the last scattering surface, and not for example the local h. 
The challenge for the WMAP analysis is first to see whether one can find at all a good fit to the CMB data, without 
Dark Energy. It turns out that one can (see section |V|, but, crucially, a reasonable fit of the WMAP data without A 
requires a relatively low Hubble parameter outside the Void: 

0.44 < /lout < 0.47, (11) 

(at 95% C.L.). 



The discrepancy between the LTB and EdS model goes like a Rees-Sciama effect, [L/th)^ , which is ~ 0(10"^), according to 34]. Such 
a correction is irrelevant for supernovae, while it could be relevant for CMB. We note, however, that [S^. |25f find a larger correction in 
the luminosity distance. The reason for the discrepancy, however, is still unclear to us. 
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Now, these two constraints (/lout from CMB and the J from Supernovae) can be combined together. And the third 
challenge now is whether we get a local value h which is consistent with local measurements of the Hubble parameter. 
Combining the range Eq. pT|) with the constraints fr-om SN Eq. (fTU|) . we get a reasonable range of 

0.51 < < 0.59. (12) 

(see fig. [7]) and we have to compare this with the local measurements. 

These local values typically vary over a wide range. The Hubble parameter measured using supernovae (35j | reads 
h = 0.59j:;^4, the Hubble Key Project [33| measures a value of h — 0.72^ Qg (although in [37| a lower value of 



h — 0.62j]Q Qg is given, with an improved treatment of Cepheids). Measurements of clusters using Sunyaev-Zeldovich 
distances [38j (which is based on data at different redshifts, up to z ~ 1) gives a much lower estimate, h = 0.541'q3 
(in EdS), as does measurement at hig h redshift (0.3 < z < 0.7) using gravitational lensing [11]: h = 0.48t;o3 (for a 
more comprehensive summary see [3). In fact, the value of h estimated also seem to decrease as one looks at sources 
with larger redshifts which would be a prediction for the MV model. However, a detailed study of this issue is well 
beyond the scope of our paper, but we want to emphasize that the local value of the Hubble parameter has a large 
window, Eq.([T]) being perhaps a fair range to consider. 

Clearly there is an overlap between Eg. p^ and Eq.(IT]), which is now consistent with supernovae, WMAP and local 
measurements of Hubble. 

This can now be used to pinpoint the underdensity contrast required in the void. As we will analytically show in 
the next section (and verify numerically), the jump parameter in LTB models does not depend on the details of the 
curvature (density) profile, but only on the amplitude fcmax, or equivalently the maximal underdensity contrast at the 
center of the void. We find that a central underdensity between 44% and 58% reproduces the relevant range Eq. pU]) 
of the jump parameter, and it is easy to check that this is also consistent with Eq.([T]) and Eq. pTj) . Notice however 
that the average underdensity is always somewhat smaller than the central value, see e.g. fig. lIVBl 

At this point one may be concerned about the plausibility of the MV model on two different accounts. Firstly, even 
if we take the observational evidence of the existence of a large underdense region seriously ^S], the underdensity 
contrast required to be consistent with WMAP and supernovae is quite large. Secondly, the local value of the Hubble 
parameter is certainly on the lower side. Both of these problems can become milder if one could obtain acceptable fits 
to WMAP with slightly higher /lout- In section |V] we briefly discuss how it may be possible to evade these problems, 
but a more detailed investigation of these issues is out of the scope of this paper. 



IV. SUPERNOVAE FITS 



A. Analytical Results 



Our aim in this section is to quantitatively fit the supernova data (we use here the dataset from ^]) using the MV 
model. In order to have better control, we decided to perform both numerical and analytical analysis. As explained 
in (34] . as long as L <C Rh, one can find excellent approximations to the luminosity distance- redshift relation. This, 
not only helps us physically understand the effects of corrections coming from inhomogeneities better, but also provides 
us with a non-trivial check on the numerical calculations. In the appendix we have obtained expressions for DL{r) 
and z{r) (which can be used to obtain Dl(z) implicitly) for any general profile. We also provide the reader with a 
summary of all the equations necessary to reproduce the analytic approximation for Dl(z) in Appendix IIX Fi in a 
self-contained form. Inside the LTB patch, the redshift as a function of the radial coordinate looks like 

z«^[l + 2/(3%)/7r)] . (13) 

while the angular distance is simply given by 

DA^r[l + f{3kir)/Tr)]. (14) 

In deriving these formulas we have used a specific choice of the radial coordinate, given in Eq.dSSl) of appendix HXB 
such that r approximately corresponds to the proper distance today. 

The luminosity distance, in General Relativity, is always related to the angular diameter distance [41] Da via 

DL^il + z)^DA, (15) 
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FIG. 1: The for Supernovae I A as a function of the jump h/hout = Hi„/Hout, for different values of the size of the 
inhomogeneous region (whose boundary ends at redshift zjump)- We have used here a model with two FLRW regions (empty 
inside and EdS outside), with two different Hubble parameters. From bottom to top the sohd curves correspond to zjump ~ 
[0.09,0.08,0.07,0.06,0.05]. The two dashed lines correspond to a 10% and a 1% goodness-of-fit. The number of d.o.f. is 181 
(we have used the Riess Gold dataset [29l]V 

and thus we now have all the ingredients to obtain Dl{z) inside the patch. One can easily verify that, in the above 
expressions for Da and z, the terms outside the brackets correspond to the FLRW results for a flat universe. / is 
an universal function (it does not depend on the profile) defined in the appendix, which gives us the deviation of 
the Dl{z) curve from the FLRW result. As one can see, our analytical results agree very well with the numerical 
solutions, see figfTOl 

Now, one defines the Hubble parameter as the initial (z=0) slope in the Dl — z plot: using this definition one can 
obtain (see appendix IIX DI for details) an exact relation between the jump parameter and the central density contrast: 

J^j^ = 2-il-\5o\)'/'. (16) 

lout 

Surprisingly, this expression does not depend on the specific form of the profile, and therefore lends generality to the 
analysis. 



B. Numerical Analysis 

We employ in this section a two steps strategy. First, without even using the LTB metric, we try to fit the data 
with a crude approximation of the void, which consists of an empty (curvature dominated) FLRW Hubble diagram 
for the inner region and then an EdS Hubble Diagram for the outer region. Between the two regions (z < zjump and 
z > ^jump) there is a discontinuous jump in the Hubble parameter i?in/i?out- In this way we get a good idea about 
what are the best values for and Zjump- The results are shown in fig. [T] 

As one can see from the plot, the larger is the value for zjump the better is the fit. However, we do not gain much by 
taking zjump larger than, say, 0.08 (which corresponds to a radius of 250 Mpc//i). It is also interesting to note that a 
2;jump as low as 0.05 (which corresponds to a radius of 150Mpc//i) still gives a reasonable fit (goodness-of-fit is higher 
than a few %). Almost independent of Zjump, the best value for the jump is around ~ 1.2. 

As a second step, then, we try to reproduce these results with a full LTB study. For simplicity we focus on only one 
value of L for the LTB patch {zjump ~ 0.085). A further observational motivation for considering such a redshift 
comes from the fact that it also approximately coincides with the redshift of the Sloan Great Wall, which spans 
hundreds of Mpc across and it could be suggestive of being the "compensating structure" expected at the boundary 
of the LTB patch Q. In the profile Eq.([5]), we therefore fix the radius L, and let fcmax vary (which corresponds to 
varying the jump J', or equivalently the central density contrast So)- 

We solve numerically for the Dl — z curve for several values of k^a.^ (which correspond to several values of J'), and 
we compute the x^- We show in fig. ^ the as a function of the jump, interpolating between the results of the 
numerics. This interpolating function is then used to compute the statistics: We find that the Icr range of the jump 
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FIG. 2: The for Supernovae lA as a function of the jump h/hout ~ H^n/ Hout, in a full specific LTB model, matched to 

_ 2 In 

FLRW at redshift zjump = 0.085. The dashed lines correspond to the Icr, 2(t, 3(t and 4(7 where we used as a likelihood e ' . 
The number of d.o.f. is 181 (we have used the Riess Gold dataset [29l]). 



corresponds to 1.214I gj^g. For the density contrast at the center this translates to 5q = 0.514_'Q3g. 

Let us comment briefly on the values that we get for the as compared to other models. The EdS model has a very 
bad fit to the data, since its for the same dataset is about 284. This has a very low goodness of fit. On the other 
hand the ACDM model has a much lower than our model ([29| quotes 150)^, which is indeed strangely too low^. 
Now, in terms of goodness-of-fit our is what one expects typically, since it is roughly equal to the number of d.o.f., 
and this makes our model a good fit to the data. On the other hand if one allows a new free parameter (Oa) then 
the best fit turns out to be at a nonzero value for JIa, and so the parameter value A = would be formally excluded 
at several a (assuming a likelihood that goes as /^). This situation is similar to what we will encounter when we 
perform the CMB fits (see section lv| : the MV model has a worse as compared to ACDM, but the question that 
we want to ask is about consistency of SN data with a MV model, and for this question the answer seems to be yes, 
the x^/d.o.f. being roughly equal to 1. 

We also note that we use only one dataset (while there are other ones in the literature), since we would qualitatively 
get very similar answer and it is not our purpose here to compare dataset with others, but just to check the consistency 
of the model. 

Finally we show, as an illustration, one example of a plot of Dl — z in figure IIVBI together with the shape of the 
density profile (as a function of z) . 



V. MCMC FIT OF THE WMAP DATA 



In order for the MV model to be viable at all, it is crucial for this picture to be in agreement with observations of 
the CMB spectra, among other things. It is commonly assumed that the ACDM model, with a non-zero cosmological 
constant, is the only one which can adequately explain the CMB spectrum. This is based on the result that once 
one assumes a "flat" prior on A, it turns out that the "most likely" parameters, given the WMAP data, correspond 
to 17a ^ 0.7. The question that we want to ask however, is about consistency of WMAP with EdS: can we get a 
reasonable fit to the CMB spectrum even after setting A to zero? To put it differently, if we had a strong theoretical 
prejudice against having a non-zero cosmological constant, or if there were other observations disfavoring it, then 
would the 3-yr WMAP data independently rule out an VI m = 1, EdS universe? (Here M means total matter ~ 



^ The open empty Universe has also a low x^> of about 160. 

® We note here that all the SN fits are plagued by not knowing exactly what are the errors on the SN measurements. In fact, if one used 
only instrumental errors, then the data points would have a very large scatter with tiny errors, and there is no smooth curve which can 
give a fit to the data. Then what is done by SN collaborations is to artificially add by hand an error bar of about 0.15 magnitudes in 
quadrature to all data points, which is typically justified saying that this is the typical variability of the intrinsic SN luminosity. This 
is what makes the concordance ACDM so low. 
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FIG. 3: In the upper plot we show a fit of the Supernovae data (Riess et al. [2y]) with an LTB model which has = 186 (the 
d.o.f. are 181). The inhomogeneous patch extends up to z ~ 0.085 and the underdensity in the center is Jcentre ~ —0.48. 
We have shown Am = m — mempty'- the magnitude (m = SLogioDi,) minus the magnitude of an empty open FLRW Universe 
as a function of the redshift z. The blue solid line is our inhomogeneous model, the red dashed- line is an EdS model (whose 
Hubble constant is normalized through the nearby supernovae), the green dotted line is the best-fit A.CDM. In the lower plot 
we show the density contrast for the same model, as a function of z. The average contrast (\/ (5^)) in the inhomogeneous patch 
is 0.43 (\/ {5^) — 0.33 in the underdensity, (S^) ~ 0.48 in the overdensity) . 



baryons -I- dark matter.) In particular, what if we introduce additional features in the primordial spectrum, rather 
than tampering with the composition of the universe? If we indeed obtain a reasonable fit using such a different 
"prior", the next important step would be to check whether this parameter set is consistent with the supernovae fit. 
This is what we plan to do in this section. 

Rigorously speaking, this question seems technically challenging because one would have to compute the secondary 
effects, i.e. what the spectrum of the CMB radiation would look like after passing through the local underdense 
region, and maybe many other such regions^, it encounters on its journey to us. According to (33 |. the corrections 
to the redshifts of photons which pass through a void of size L is a Rees-Sciama effect that goes like {L/RhY. A 



The assumption that we live in a void could naturally lead us to consider that the universe might contain many such voids, a bubbly 
universe. In this case one would have to compute the passage of the photons through many such voids. 
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TABLE I: Priors for different parameters in the COSMOMC Runs. Here flth^M is the physical baryon density, fim/iout is the 
physical dark matter density, Zre is the redshift at re-ionization, Us is the spectral index, Os is the running of the spectral index 
and As is the amplitude of scalar fluctuations (for deflnitions see, e.g. [l^). 



coherent addition of this effect due to many voids could produce a correction of order {L/RhY ■ Thus for a void with 
a typical radius ~ 200/ft- Mpc that we considered in this paper, such a cumulative effect could be ~ 10^^ — 10^"^. 
This can be ignored for the study of supernovae*. On the other hand, if these many voids exist, they would give a 
sizable effect on the CMB. The number ~ 10~^ — 10"'^ would refer to a monopole in the CMB, while the correction 
to higher multipoles would be smaller (depending on how different is the number of voids along different directions 
in the sky). However, in this paper we ignore such secondary effects. On the qualitative side, in fact, we expect 
this to be important only for small I and decay fast for larger I, and it should act in the same way as an Integrated 
Sachs- Wolfe effect ^. 

The correction to the CMB redshift that comes from our local void, instead, will depend on how symmetric the void 
is, and how "centrally" we are located. For an off-center observer, in appendix IIX El we perform a non-perturbative 
estimate of the dipole moment, and find that in order for it to not exceed the observed value ~ C(10~^), "we" must 
be located very close to the center, approximately within 10% of the void-radius, in concordance with the findings 
in [i^. In this case the correction to the higher multipoles are much more suppressed and not visible in CMB [43 |. 
Departure from spherical symmetry, on the other hand, may have a much more interesting effect, specially on the 
lowest Is in the CMB spectrum, and could be visible^°. However, such a study is clearly out of the scope of the present 
paper; instead we will restrict ourselves to spherically symmetric voids and neglect these possible secondary effects 
on the CMB coming from the voids embedded in the homogeneous EdS background. Thus, the question reduces to 
whether the CMB spectrum can be reproduced given an EdS background. 

As one would expect, we find that if one assumes as priors, no dark energy, as well as no additional features in the 
primordial spectrum (other than spectral index and amplitude) , one obtains a very poor fit to the 3-yr WMAP data 
(see table IIV|) . However the situation changes if we introduce a possible "running in the spectral tilt" , as , in the 
observable ~ 7 e- folds of our universe in CMB (following the same definition as in [iy|) 

We have performed a Monte Carlo Markov Chain (MCMC) analysis of the WMAP 3 year data using the program 
COSMOMC [IB]. Our runs were performed with the priors given in table fVl We used the version of the COSMOMC 
program which lets one analyze the range 2 < ^ < 30 for TT correlations and the range 2 < Z < 23 for TE-I-EE 
correlations using the pixel-based approach (T, Q and U maps), which offers a much more accurate treatment of the 
low-/ likelihood [43|. One has (957-1-1172) pixel data in all. The rest of the correlations that we considered consisted 
of Cf^ in the range 31 < ? < 1000, and Cf^ in the range 24 < Z < 450. 

We find that an EdS universe with no dark energy but with a value of the Hubble constant, -ffout, significantly lower 
than the conventionally accepted value of 70 km/s/Mpc gives a very reasonable fit to the CMB spectrum, see fig. HI 



We mention, again, that [24l |25|| find a larger correction to the luminosity distance, which could be potentially important for super- 
novae [24j . Since the reason of the discrepancy is still unclear to us, we do not discuss it here. 
This might explain the claimed detections of ISW correlations [40| . even without invoking Dark Energy. 

In this context we note that similar effects in anisotropic geometric void configurations have been used to explain the low multipole 
anomalies in the CMB sky Il2l . 

Ideally, one would like to introduce two additional scales where significant running of the spectral tilt starts and ends. This would also 
obviously improve the fit. However, to keep the analysis simple, we have assumed that these two scales lie outside the observed spectrum 
in WMAP. 
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FIG. 4: ACDM and EdS fits to the WMAP 3 binned data 
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TABLE 11: Most likely parameter values with 1 a errors for the various COSMOMC Runs 



For the high multipoles (31 < I < 1000) TT power spectrum our goodness-of-fit (G.F.) is around 2%, compared to 
around 5% of the concordant ACDM model. For the overall fit of both the TT+TE+EE spectrum involving 3520 
d.o.f., the EdS model has a reduced^^ Xeff r of 1-016 with a 26% G.F., as compared to the "concordance" ACDM 
model^'^ with xltf r — 1-005 and G.F-=41% (see table llVl for more details)- The most likely parameter set along with 



The "effective" is obtained directly from COSMOMC 16]. To obtain tlic reduced effective cfii-square, Xett r' ™^ J^^* divide it -with 
the number of independent degrees of freedom. 

The "concordance" best fit ACDM model is obtained by running the COSMOMC program including both the WMAP 3-yr and 
supernovae data. The best-fit ACDM parameters for the WMAP 3-yr data alone yield very similar Xett n indicating the presence of 
the well-known degeneracy in Ojvf ~ ^ plane of the WMAP data. In fact, it is this degeneracy that we exploit to fit CMB with Qm = 1 
and low Hubble parameter. 
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FIG. 5: Parameter likelihoods to the WMAP 3-yr data for the run "EdS, as 7^ 0". Dotted lines are "mean likelihoods" of 
samples, while solid lines are "marginalized probabilities" p^ . 

their Icr bounds are tabulated in table HIl also see the likelihood plots, fig. [H 

We also produce two 2-dimensional likelihood contour plots: (i) /lout vs. fim/flf, which are the only two independent 
parameters related to the composition of the universe, and (ii) rig vs. as which characterize the spectrum. 

The most crucial quantity to consider is the Hubble parameter and in particular what a consistency with the supernova 
data implies for the locally measured value. In fig. [7] we show a contour plot combining the constraint from supernova 
fit in the previous section with that of WMAP. As promised before, we find that the locally measured Hubble parameter 
can be as high as /i ~ .59 at the 2<t, or 95% C.L., which is within the acceptable range of the different measurements 
of the Hubble parameter. 

Let us briefly discuss about the values that we obtain for the other cosmological parameters, a more detailed discussion 
on some of these constraints is presented in the next section. The main constraint on the baryon density comes from 
BBN, and we are indeed consistent with the data (see next section for details). As one can see from the likelihood 
plot, fig. [5] as well as table [TTl the ratio between dark matter and baryons is somewhat higher, Um/^b ^ 10, than the 
"concordance" ACDM model value of flm/^b ~ 6. Measurements of light-to-mass functions in galaxy clusters can 
in principle be used to constrain these numbers, but presently they suffer from relatively large uncertainties (see for 
instance [44| | , and references therein) . The issue is further compounded by the fact that our local ratio of abundances 
within the LTB patch may not represent the global ratio. A more detailed investigation will be required to settle the 
issue, but potentially this could be a problem. For the total matter density, one now has tight constraints from the 
observation of BAO ^45]. As we discuss in the next section, the total matter density in our model (which is the same 
as the critical density and hence cx h^^^.) seems consistent with these measurements. 

What about the properties of the primordial spectrum? Our best fit spectral tilt is relatively low, rig ~ .73, but 
there are several inflationary scenarios where such low spectral tilts are common (for example in modifications of the 
old inflationary scenario from false vacuum [48| , or inflation from exponential potentials naturally occurring in string 
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FIG. 6: Contour marginalized likelihood plots to the WMAP 3-yr data for the run "EdS, Os ^ 0". The coloured map 
corresponds to mean likelihood, while the solid lines correspond to marginalized 1-a and 2-a contours. 
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FIG. 7: 1-a and 2-a Contour plots for h vs. hout- The blue bands come from the SN-I analysis, while the red bands correspond 
to constraints coming from WMAP. 
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TABLE III: Best-fit Minimal Void Model Parameters 



theories, see for instance [49j). Our model also requires a significant running, ag ~ —.16. It is a known fact that 
the 3rd year WMAP data favors a significant running of the spectral index which deviates from a Harrison- Zeldovich 
scale invariant scalar power spectrum. For example, the analysis of 50] gives a running as ~ ^^■^^^t.^o'ff^ 60% 
confidence level. In fact, most inflationary models predict a running spectral index 51,1 (see also ^52|]; models of 
inflation from a False Vacuum have typically an abrupt transition in the spectral index [48|) . Additional constraints 
on {ns,asT<J%\ can mostly come from observations of large scale structure and weak-lensing experiments. In the 
context of our MV model, this is a difficult and somewhat tricky task which we have postponed to a future analysis, 
however we do discuss briefly possible implications in the next section. 

Finally, we note that our value of the re-ionization epoch (optical depth) is broadly consistent with the usual obser- 
vations [5^1 (see also discussion in [1^). 

To summarize, our best fit (WMAP -I- SNIa) MV model consists of 8 parameters, one of which, the length scale of the 
void, has been chosen at the value L — 250 /h to derive our best-fit model. However, as noted in the introduction, if 
one "accepts" a G.F. ~ 10% to the supernovae data, then one can go down to a much smaller length scale, L ~ 160/ h. 
Out of the other seven parameters, six of them (columns 2 to 7 in the Table of IIIII are obtained from the fit to 
the WMAP 3-yr data using COSMOMC, while the last one, (column 8), is constrained from the supernovae data. 
We note that a "minimally acceptable" model with respect to the central underdensity contrast would be obtained 
with a maximally acceptable /lout ~ 0.47, at the 95% C.L.. This in conjunction with Eq.((T|), then tells us that the 
minimal jump parameter has to be 1.17, or equivalently Sq ~ 0.44. Using these information we tabulate all the 
parameters in Table [Till for our "best-fit" and "minimally-acceptable" model. We note that the values of Sq and L in 
the "minimally-acceptable" fit is not far from what observationally is suggested in . 



VI. CAN WE IMPROVE WMAP AND SUPERNOVAE FITS? 

We have seen that by allowing significant running in the range of the observed CMB spectrum one is able to obtain a 
reasonable fit to the WMAP 3yr data. However, the overall fit is not as good as the best-fit ACDM model. Secondly, 
as is clear from the combined contour plot fig. [3 consistency with WMAP and supernovae data requires a relatively 
low local value of the Hubble parameter. The underdensity contrast required is also quite high (centrally around 50%, 
and on average around 35% in the Void). Can we somehow modify the MV model to get a better fit and overcome 
these difficulties? We now discuss two different modifications in this context. 



A. "Bump" Model 

The first one concerns using different "priors" for the primordial spectrum. For instance, in p7| the authors assumed 
the existence of a bump in the primordial spectrum as a prior, rather than considering an overall running as we do, 
in order to fit the CMB data without Dark Energy. Although in these models the number of parameters is larger 
than what we consider, one obtains much better fits to the WMAP data (in fact, slightly better than ACDM), and 
is thus worth investigating further. Such a bump can be produced by a rapid succession of two phase transitions (l7j 
and is thus phenomenologically well-motivated. Moreover, it is rather intriguing and promising to note that [l^ such 
a bump would also enhance the probability of having voids today at the scale of the bump itself, which happens to 
be approximately the same scale we are considering here. This "bump" model, in its original form, of course cannot 
reproduce the supernovae data, and the Hubble parameter (/lout ~ 0.44) is too low. So it seems natural to merge this 
model with our MV scenario. Can the parameter set obtained be consistent with the supernovae analysis that we 
have performed using the local void? 

Of course, having a local void again ensures that the supernovae data is consistent. The crucial question is whether 
putting together the MV framework with the "bump" model could lead to an "acceptable" local Hubble parameter. 
As we see in the contour plot (see fig. [5]) at the 95% C.L. one can have as high as h .57, which is definitely within 
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FIG. 8: 1-a and 2-a Contour plots for h vs. hout for CDHM bump model [T3|. The blue bands come from the SN-I analysis, 
while the red bands correspond to constraints coming from WMAP. 

the acceptable range Eq.([T]). 



B. Adding Curvature 

In this subsection we consider a different possibility, namely adding curvature to the model. Although having curvature 
would be considered fine-tuning in an inflationary paradigm, we point out that the low multipole anomalies [55^, if 
taken seriously, could be suggestive of having only the "minimal" number (50-60 depending upon the reheating 
temperature) of efoldings, which would be consistent with having a slight curvature. Also, we note that other models 
involving cyclic scenarios typically do not predict a flat universe to any high precision. 

Accordingly, we performed a run where we allowed up to 5% in curvature along with including running of the tilt, 
as before. We found that the best-fit parameter set prefers the highest value of spatial curvature that we allowed. 
Consequently, we performed a run with flk = 0.05, corresponding to a slightly closed universe to see how curvature may 
affect the goodness of fit-^"*. We now indeed find a much better fit to the WMAP data. For the overall TT-I-TE-I-EE 
data, Xeff,r = 1-012 corresponding to a 31% goodness-of-fit (see table |IV] for more details). The Hubble parameter, 
however is slightly lower than our previous results, as can be seen from the likelihood plots involving (ii) Qk and /lout, 
also see table [TTl None of these results are very surprising or new. Previous studies had already observed that one 
can get good fits to WMAP with a closed universe, but it is precisely because of the rather low value of the Hubble 
parameter required for these fits that these models are not considered seriously. However, when combined with the 
jump parameter obtained from the supernovae analysis^^ given by Eg.ljlOp. this gives us a local Hubble parameter 
which can be consistent with observations, given in Eq.((l]). 

It is also worth pointing out that, as is clear from the contour plot fig.O there is a degeneracy direction in the WMAP 
data where as we simultaneously increase the curvature and the Hubble parameter we can still get good fits. This 
suggests that even if we allow for a slightly closed universe, by decreasing the Hubble parameter slightly (from our 
EdS value) we may be able to get significantly better fits to the WMAP data. In other words, when combined with 
other data, such as measurements of local Hubble parameter which prefer higher values of the Hubble parameter, the 
MV model may still provide a reasonable fit. 

In passing we note that, the 2a range for the tilt and the running is much closer to the conventional values as compared 



We are currently pursuing a more exhaustive analysis of the void model with curvature. 

In principle once one adds curvature, one has to redo the analysis of the supernovae data set. We have not done it for this preliminary 
analysis, because we do not expect any significant difference from the small amount of curvature that we allow. 
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FIG. 9: Marginalized Likelihood plots for the WMAP 3-yr data for the run "EdS, as, 7^ 0". The coloured map corresponds 
to mean likelihood, while the solid lines correspond to marginalized 1-cr and 2-(j contours. 
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4.7% 
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0% 
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3652.3 


6% 


EdS / 
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1.9 % 


1475.5 


5.7% 


3577.4 


24.6% 


EdS Qs, fife / 


1048.7 


2.9% 


1466 


7.9% 


3560.9 


31.1% 



TABLE IV: xiff and goodness-of-fit for the different COSMOMC Runs. The first column corresponds to high-/ TT power 
spectrum, (31 < / < 1000). The second column corresponds to both the high-/ TT (31 < / < 1000) and TE (24 < / < 450) 
data. Finally, the last column contains the total statistics of TT {2 < I < 1000) and TE (2 < / < 450) spectrum. 



to our original MV model (see figs. Eland [H for comparison). 



VII. CONSISTENCY WITH OTHER OBSERVATIONS 



Here we briefly discuss the consistency of our MV model with observations other than the supernovae, WMAP and 
local measurements of the Hubble parameter. 



BBN: Primordial Nucleosynthesis has been a spectacular success story for the Standard Big Bang paradigm which 
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predicts specific freeze out abundances of light elements such as D, He^, He^ and Li. These freeze out abundances 
depend on the baryon-to-photon ratio. Since from the measurements of CMB temperature we know the photon 
energy density precisely, BBN can also constrain the baryon density in our universe today. The success of the BEN 
paradigm lies in the general agreement of this number, measuring the abundances of the different light elements 
spanning 9 orders of magnitude (for a review see [54|). BBN therefore constraints the baryon density, so that at 95% 
C.L. we should have-'^^ 0.017 < 17^/1^^^ < 0.024. This is indeed consistent with the parameter range that we obtain 
from the WMAP run, ^^K^^,^ = O.OISI qqj. It is remarkable that although we have a higher baryonic abundance, the 
lower Hubble parameter almost precisely compensates to yield approximately the same baryonic energy density as it 
is obtained in the "concordance" ACDM model. 

BAO: Recently, a remarkable achievement of observational cosmology has been to identify the baryon acoustic peak 
in the galaxy-galaxy correlation function using Luminous Red Galaxies (LRG's) [1^. The overall shape of the galaxy 
correlation function mainly depends on the shape of the primordial spectrum (tilt and running) and the epoch of 
matter-radiation equality (scales which entered the Hubble horizon before the equality have their amplitudes relatively 
suppressed as compared to the ones which entered later). On top of the "overall envelope" one has now observed 
a tiny peak coming from the baryon acoustic oscillations. The position of the peak is related to the sound horizon 
of the baryon-photon plasma at the time of recombination. In fact, what one really measures is more like an angle 
which is the ratio of the sound horizon at recombination (evolved at z « 0.35, which is the average redshift of the 
LRG survey) and the angular distance^^ at the same redshift z ~ 0.35. This ratio therefore is not only sensitive to 
the baryon density in the universe, but also to the evolution of our late-time universe and therefore, to the amount 
of dark energy, for instance. Using essentially the two pieces of information (overall shape and peak) one is able to 
constrain two different quantitieSjfor instance, the matter density and dy (a specific combination of the transverse 
and angular distance at z « 0.35 '45'|). This in turn can constrain the composition of the universe and it was claimed 
in [4^ that a pure EdS model is ruled out at the level of 5(7. Can the MV model be consistent? 

Firstly, it is difficult to provide a crisp answer to that question based on the analysis done in [45| because the analysis 
of the data (conversion from redshift to distance etc.) is done using the "concordance" ACDM model. In particular 
we point out that precisely in the redshift range of the sample, 0.16 < z < 0.47, the luminosity distance vs. redshift 
curve of the void model (which is the same as an EdS model, in this range) differs significantly from the ACDM curve. 
Thus to be precise, one needs to reanalyze the LRG data in the context of an EdS model. Nevertheless, one can try 
to see whether one can satisfy the bounds on dv and Qmh^ that was placed in jisj : 

dy = 1370 ±128 and fi^/i^ = 0.130 (n^/0.98)^-2 ± 0.022 , (17) 

where the errors correspond to approximate 2<t (95% C.L.) values. Now, in our model, we have a low spectral tilt 
(and also a relatively large running, which can alter the shape of the correlation function and hence the constraints). 
Relegating a more systematic analysis for future, and just correcting for the lower tilt in our model implies the 
following constraint for the total matter density, which is given by the average Hubble parameter in our model: 

rj^/i^ ^ /i^^t = 0.185 ± 0.022, (18) 

where we have used our best-fit spectral tilt, Ug ~ 0.73. We recall that our best fit Hubble parameter gives h^^^ ~ 0.205, 
and therefore is consistent with the above bound. 

On the other hand the angular distance at z = 0.35 for our model does not appear to be consistent with the values 
reported in [45|| . In fact one can check that an EdS model has roughly the same distance of a concordance ACDM model 
at 2; = 0.35 if the ratio of the value of the Hubble constants of the two models is around 1.2. Since the concordance 
model (which fits the BAO scale) has h « 0.7, an EdS that fits this scale should have hout ~ 0.7/1.2 « 0.58. As 
we have discussed, this value is too large with respect of our analysis of the WMAP data. More work is needed in 
order to find whether it would be possible to overcome this potential problem: for example adding more curvature 
could give a higher /lout, from WMAP (which, by the way, would make the whole scenario in better agreement with 
other data as well: for example with the local measurements of h). It has to be seen, through a combined statistical 
analysis including the BAO data, whether this could give a consistent picture. A different (though not very appealing) 
possibility, which would certainly work, is to make the Void much larger, extending up to redshifts of order z « 0.4. 



While estimates from He'^ and D are slightly higher, measurements using Li suggests a lower number. 

In the survey, one really measures a combination of the angular and transverse distances, see [4^ for details. 



18 



We have checked that this can give the correct distance at z « 0.35 (see also the recent analysis in fi^), and that 
moreover this gives also a good fit of the BAO scale at z w 0.2, given in 

Observations from Large scale Structure and Weak lensing: An important class of cosmological observations 
comes from large scale structures and weak lensing observations. These typically produce constraints on as, as well 
as on the shape of the primordial spectrum, Us and as (for instance using Lyman-a forest). However, as mentioned 
above, one has to revisit these analysis in the light of MV model, as one has a non-standard Dl{z) relation. We leave 
for future such a careful study of the large scale structure, Lyman-a and weak lensing data. Let us still make a few 
brief comments. 

About (Tg, at first sight our value is a bit high, erg ~ 0.92^ q8j model, but even if this situation turns out to 

be incompatible with the large-scale structure data (after a careful study) , this may only be indicative of the need to 
include some hot dark matter component [TtI . [isj . In the light of neutrinos having mass, this is a perfectly natural 
scenario to consider. About Ug and \as\, the values are respectively lower and higher than what the conventional 
ACDM analysis suggests and in particular one may worry about confiicts with Lyman-a measurements. However, 
we firstly point out that an analysis of the Lyman-a measurements has to be now re-done in combination with the 
different set of priors that we use to study the WMAP data. Secondly, introducing new physics, such as including a 
little curvature, can push the values of and a^ much closer to the standard values. In short, there are too many 
uncertainties for us to make here any concrete conclusions, and one really needs to perform a careful study of the 
above mentioned observations. 

ISW Correlations: Another interesting piece of evidence for Dark energy is given by the Integrated Sachs Wolfe 
effect, which is claimed to be detected with some significance by some collaborations [i^. The detection is a correlation 
between the CMB maps of the sky and the galaxy surveys, which cannot be explained in an EdS universe (since in this 
case the linear gravitational potential does not evolve and therefore CMB photons do not get any net frequency shift 
when passing through a potential well), and therefore are interpreted as independent evidence for Dark Energy (since 
the potentials can evolve in ACDM). However the effect is absent only at the linear level, and it exists also in EdS 
in the presence of nonlinear gravitational clustering. This is usually assumed to be smaller then 10~^, but it actually 
happens to be of order 10^^ (and thus, visible in the CMB) for structures as large as those that we are proposing in 
the present paper (few hundreds of Mpc) . It would be interesting to try and reproduce the ISW detection assuming 
the presence of large voids and structures in the sky. 

Moreover in the local underdense region we have assumed that the growth of fluctuations is different than the flat 
CDM model (it is in fact more similar to an open Universe): this leads also to an ISW effect for density fluctuations 
localized inside the Void. Studying this effect would be very interesting and could significantly affect the \ow-l part of 
the CMB spectrum and therefore also the parameter estimation from the CMB. However this goes beyond the scope 
of the present paper since it would require a full treatment of the growth of density fluctuations in an LTB metric 
(this problem has been recently attacked by [5!|). 



VIII. CONCLUSION AND DISCUSSION 

The Type la supernovae data reveal that our universe is accelerating at redshifts that approximately correspond to 
the epoch of non- linear structure formation on large scales (the epoch of the formation of the so-called "cosmic web" ) . 
Given this fact, we have explored the possibility that the effect of a large scale void can account for this acceleration 
due to a jump between the local and the average Hubble parameter, instead of invoking a spatially constant dark 
energy /cosmological constant component. We find that the Minimal Void (MV) model can consistently account for 
the combination of the Type la supernovae, WMAP 3rd year, BBN constraints, provided that the void spans a radius 
of about of 200 Mpc/ft. with a relative under density of 45%, near the center. The MV model can accommodate 
reasonably all of the data considered, although the fits are not as good as the concordance model. However, we see 
the possibility of obtaining just as good fits when one includes curvature or invokes non-standard features in the 
primordial spectrum (a "bump" for example). We leave these issues for an upcoming work. On the other hand we 
have seen that the Minimal Void is in trouble with the Baryon Acoustic Oscillations measurements, since outside the 
Void, the Dl{z) curve is just the usual EdS one, and the Hubble parameter /lout from WMAP is too low. More work 
is needed in order to find whether it would be possible to overcome this potential problem (for example by finding a 
fit for WMAP with higher hont)- 

We end with observational and theoretical possibilities of distinguishing the MV model from ACDM. The MV 
model predicts that the spectral index has to run significantly in the WMAP3 data and that the "average" Hubble 
constant {i.e. outside the local region) has to be around /lout ^ 0.45. The AC DAI model, instead, requires a finely 
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tuned cosmological constant or dark energy component, in order to be consistent with the same data set. Both cases 
require significant model building and new physics that are currently being pursued by the community. How are 
we to distinguish between these two models? The first logical way seems to perform galaxy counts up to very large 
distances and in a wide area in the sky, in order to directly check if we could really live inside a huge Void. Moreover, 
there are features which can be checked by looking at SN la themselves: firstly, the luminosity-redshift curve in the 
two models deviate from each other significantly at redshifts z > 1. Secondly, in the MV model the curve has a 
sharp peak (in correspondence with the boundary of the local region) around z ~ 0.1, while this peak does not exist 
in the KCDM model. The up-coming experiment SDSS-II [S^] will probably be able to discriminate the presence 
of such a peak. Another unique prediction for the MV model comes from realizing that the void is not expected to 
be exactly spherically symmetric, which could lead to detectable anisotropies in the Hubble parameter as well as in 
the low multipoles in CMB. Additionally, these anisotropies should be correlated! We note, also, that one could be 
able to constrain Voids by looking at the blackbody nature of the CMB [H, [s^l . Our MV is still consistent with 
these constraints (while, according to [s^l, voids that extend up to z ~ 1 are excluded). Finally, studying large scale 
structure (as we plan to do in future work) one can study the compatibility of the primordial power spectrum we 
are assuming (with low tilt and large running, or with a bump) with the matter power spectrum. It may also be 
possible to test the existence of such a large running using Planck-satellite data as suggested by the Bayesian analysis 
performed in [sl] using simulations. 

In conclusion, we have shown that, for WMAP and SNIa observations, the MV model could be taken as an alternative 
to invoking a dark energy component that will be further tested in forthcoming supernovae observations. On the 
other hand this has to be made consistent also with the Baryon Acoustic Oscillations. On the theoretical end, much 
work needs to be done to establish if such large voids can actually be produced in our Universe by generic physics of 
structure formation. We are currently pursing this issues. 



Note Added: Most of the above research work was completed before the release of the WMAP Syr data 
and we have decided not to re-analyze the CMB data in the present paper for the following reason: although the 
5-yr data improves the 3-yr data, there is no significant qualitative difference between the results presented in the 
3-yr and 5-yr survey. In this context, we further emphasize that our aim in this manuscript is not to compete with 
ACDM on the basis of Bayesian likelihood analysis (in which case the analysis can be very sensitive to the data, for 
instance a difference of ~ 1 — 2 may be significant), but to simply present a model which can be consistent with 
the data on the basis of the goodness of fit (for instance, a difference of ~ 1 — 2 does not significantly reduce 
the goodness of fit). In addition a more systematic treatment including other cosmological data (BAO, Large Scale 
Structure data) and more recent data (CMB and Supernovae) is the subject of a future publication. 
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IX. APPENDICES: ANALYTICAL RESULTS FOR LTB METRIC 
A. Metric & Density Profile 

In our paper we are interested in a special class of exact spherically symmetric solutions of Einstein's equations with 
dust, known as the "open" LTB metric (in units c—1). We follow the treatment given in ([3, 34]), where we have set 
the "mass function" to be cubic, which amounts to a redefinition of the radial coordinate (which is always possible if 
the mass function is a growing function of r) . The metric is given by: 

ds^ = -dt^ + S^{r, t)dr^ + R^{r, t){de^ + sin^ Odif^), (19) 
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Here we have employed comoving coordinates (r, 0, if) and proper time t. The functions ^^(r, t) and the dust density 
p{r, t) is given in terms of -R(r, t) via 

where a dot denotes partial differentiation with respect to t and a prime with respect to r, while the function R{r, t) 
itself is given implicitly as a function of {r,t) via an auxiliary variable u{r,t): 

R{r,t)^^^(coshu~l), (22) 

^'^^it=^^^{,in\,u-u), (23) 

In the above expressions, the "curvature" function k{r) is left arbitrary (except that k{r) > 0) and this is what 
controls the density profile inside the LTB patch, while M is just an arbitrary (unphysical) mass scale. Also, we have 
introduced the dimensionless conformal time r for later convenience. 

We also note that the average density inside the LTB patch is equal to the outside FLRW density (see for instance [1, 
[S^l), in the limit in which we can neglect (Mr)^fc(r) in Eg. ([20)1 in the spatial metric when defining the average (in 
our case the correction is always negligible) . 

To get an intuitive and analytical understanding of how the density profile is related to the curvature function it is 
instructive to look at the "small-w" approximation where we only keep next-to-leading terms in Eq. (|22p and Eq. 
This gives us Eq.©. 



B. Photon Trajectories 



In order to perform supernovae fits we need to compute the luminosity (or angular) distances and redshifts for a 
photon trajectory emanating (backwards in time) from the central observer. The first step in this direction is to solve 
for the photon trajectory: 

ds^^0^'41^-^£M)^. (24) 
dr y/l + 2(Mr)2/c(r) 

The negative sign in front takes care of the fact that the time increases as the photons go towards the center. Analytical 
progress in solving the above equation is possible by realizing two things. Firstly, all quantities (t(r), z(r), Di(r)) 
can be expressed as a power series in, Mr ^ r/Rn, and since this is a small quantity for the relevant inhomogeneous 
patches, we can just keep the next-to-leading order terms in these expansions ^34]. Secondly, formally one can combine 
Eq. iP^ and Eq. ipSll to give us a power series expansion for R{r,t) explicitly in terms of {r,t) [3^ : 

R[r, t) = ^tt^\t^ (1 + R2ul + Ri4 + ...) = ^ttj^t^ (l + /(u^)) , (25) 



where 

Mo = 7(Mi)'/^v/%yand7= ( — ) . (26) 



9V2 



It is important to realize that the coefficients {Rn}, and hence the function / are universal (do not depend on the 
specific curvature function). It is implicitly defined via 



1 + f{ul) = — 11 and 6{smhu - u) = uf^ . (27) 



21 



This is what ahows us analyze the problem in its full generality. 

It is convenient to recast the equation in terms of the conformal time, r, and the dimensionless radial coordinate 

f = Mr. (28) 

Substituting Eq.lgSl) in Eq.dH]) one finds 

dr _ f 7' [1 + Er i?2„7'"T2»(ffc")'] 



dr 



Vl + 2fcf2 



(29) 



The prime now denotes differentiation with respect to the rescaled r. This can now be solved perturbative in r to 
give us 



(30) 



The first two terms within the brackets corresponds to the FLRW expression for the trajectory while the rest of the 
terms give us the largest corrections coming from the inhomogeneities within a local patch. For corrections outside 
the patch see [s^l ■ By comparing with Eq. (P5)l the above expression can succinctly be written as 



T(r) = TF{r) - -^l^rf{j'^T^k{r)) 



(31) 



where the subscript F corresponds to FLRW. 



C. Luminosity Distance vs. Redshift 

Having found the photon trajectory, the next step is to compute the redshift which is governed by the differential 
equation [2] 



dz _ (1 + z)R' 
dr ~ Vl + 2fcf2 ■ 

Again, if we are only interested in computing corrections up to linear order in f, then the redshift is given by 

df 



^'^[l + ;^(n + l)i?„7^V"(.fc")'] = ^ 



1 + z 



dfY,{n + l)Rnl^''T^"^\rk"y 



(32) 



(33) 



To evaluate the first integral we note that we can replace r by tf as we will only be making an 0(r^) error. Thus we 
have 

df I" df ^ f '^^^ ^^"^ 1 

T j Tp 7r72 J TF TT To 

The second term can be integrated straight forwardly up to linear terms in f: 

^(n + l)i?„72" j T^^'-^rk'^ydf « ^(n + l)R„-f^'W^''-^ J (fk'^ydf = ^(n + l)i?„72"r2"-iffc" 



where we have defined 



Putting everything together we have 



dfix) 



dx 



l + z = 



TF{f) 



exp 



f 2^7^f[/(7Vo^fc(r)) + 7Vgfc(0/i(7Vo^fc(0)] 
I 9to 



(34) 



(35) 
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Thus we have obtained an analytical approximation for the redshift as a function of the radial coordinate. We note 
in passing that the term in front of the exponential precisely correspond to the FLRW result. The corrections come 
from the exponential. In fact for small z one finds 

z « |^7'^[1 + /(7Vo'fc(r)) + 7^0^^(0/1 (7V2fc(r))] . (36) 
9ro 

The luminosity distance, in General Theory of Relativity, is related to the angular diameter distance. Da via 

i^L = (1 + zfDA ■ (37) 
Now, in an LTB model when the observer is sitting at the center, the angular distance is simply given by 

Da^R^ \n^\T^ (1 + f{l^Tlk{f))) . (38) 

Thus we now have both the luminosity distance and the redshift as a function of the radial coordinate and one can 
easily plot Dl{z) and check whether the local void model can provide a good fit to the supernova data or not. 



D. The "Jump" 



A particularly important quantity that can be inferred from the Dl{z) curve is the jump parameter, defined by 
Eq.Q. Surprisingly, this turns out to not depend on the specific profiles, let us here see this analytically. First 
observe that since k' vanishes at r = 0, we have the general result 

i?'(0,i) = i7r7V2(l + /o), (39) 

where /o corresponds to the value of / at r = 0. Then using the exact expression for the density function Eq. (|2ip one 
finds 

67rig(l + foy 

The underdensity contrast at the center, Sq now can be easily related to /q: 

So = (l + /o)"' - 1 ^ 1 + /o = (1 + <5o)-^/3 . (41) 

Now, on the other hand using the definition of the Hubble parameter Eq.®, the correction to the redshift Eq. . 
and the luminosity distance Eq. ([38| one finds 

r-l rr-1 l + /0 



Hn ^ — H^„l- 



° °"M + /o + ^i2/l,o■ 

Or in other words 

j^J^ = l±A±^. (42) 

"out 1 + Jo 

Since 5q uniquely determines /o via ([1T|) . and /(ug) is a given function, it also determines Uq and /i,o = /i(uo)- Thus 
in turn it also determines the jump parameter uniquely. 



E. CMB dipole moment 



Let us consider our observer to be located slightly off-center, at r = tq. In this case the non-zero radial velocity of 
the observer will contribute towards a dipole moment in CMB: 

ST 

— ^ do , (43) 
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where the proper radial distance, do, of the observer is given by 



do = dr 



y^l + 2(Afr)2fc(r) 

Now, in our profile k{r) remains almost a constant for almost the entire underdense region. Assuming we are living 
in this "constant" underdense region, we have 



27r(coshu-l) r° dr 27r(coshu - 1) sinh ''{My/2k^ro) 

do ~ - ^ 



(The simplification occurs because u and hence R' becomes only a function of time.) Further, since Mro is expected 
to be very small, we have 

do = ''<^^'\^-^yo . (44) 

Taking the time derivative and simplifying we find 

• doH^M MgSinhu 

do — T — 7^- (45) 

4 (coshM-l)2 ^ ' 

We now note that u{uq) is a known function Eq. (|27p . in turn uq is known in terms of Sq via the function /(wq), see 
Eq. ()4ip . Thus, in principle, the second term in the right hand side of Eq.(|45p is determined in terms of the central 
underdensity contrast. Also, since the measured value of the CMB dipole moment ~ 10"'^, naturalness arguments 
suggest do to be of the same order, and thus we have (after some simplifications): 

doffout-lO-^ , x/2(l + /o)2 

For voids of around 200/h Mpc, and central underdensity contrasts between 40% and 50%, the dipole constraint 
Eq. (j46l) typically imply that "we" have to be located within 10% of the void radius. 



F. Analytic expression for the Dl — z curve 

In this subsection we wish to provide the reader a self-consistent summary of all the equations which are needed to 
plot the I?L — z curve, in an analytic form. Following this, a fit of any experimental dataset can easily be performed. 
Here is the set of equations, which give and z as a function of the radial coordinate r (therefore implicitly Dj^ — z). 
First of all one needs to define the function /(mq), implicitly given by: 

^ ^(cosh(7.) - 1) _ 
^ ~ 32/3(sinh(u) - u)2/3 > 

uo ^ 6i/^(sinh(u) -u)i/3^ (48) 
Then, one can use this function in the following equations: 

T(r) = TO - ^j^Mr[l + /(7V2fc(r))] , (49) 



1 + ^(0 = ( 7^ ) exp 



2 r.__2 



(50) 



Ddr) = ^^\T{rf[l + f{^\lk{r))][l + z{r)f (51) 



9\/2\ 



(53) 
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Zjump =0.085 ; 5cENTRE=-0.48 




FIG. 10: Comparison between analytic and numerical Dl — z curves. The numerical curve is the blue solid line, the analytic 
approximation is the black short-dashed line. We have plotted also the EdS curve (red long-dashed line) and the A CDM, 
with = 0.7 (green dotted line). We have used the value L = 400, with the units given in Eq. (|55p . and fcmax = 2.2 (which 
corresponds to a density contrast at the center Jo = —0.25). 



The above formulas are completely general for any LTB profile, but we now focus into our specific one given by 

2 



k{r) 



1 



L 



(54) 



Then one has to choose appropriate values for iJg, and for the length units for the coordinate r (given by M). A 
simple choice is to set: 



(55) 



y Af = iJout = /lout/3000 , 



where we have chosen, in this way, the units Mpc=l (which turns out to be a convenient choice for the problem). 
Once this is done the physical parameter L (the radius of the patch) is approximately given already in Mpc. The 
comparison between the obtained curve and the fully numerical curve is shown in fig. ()10|) 



Finally the reader may play with the two parameters: the size L and fc^ax (which sets the amplitude of the density 
contrast). We also recall that the density profile is given by Eq.([2]) and that fc^ax can be directly related to the density 
contrast at the center of the void at the present time, via the following equation: 



'5o = [1 + /(7%'fcmaxr' - 1 



(56) 
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